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Abstract 

We review recent progress in determining the partition function of the ABJM theory 
in the large N expansion, including all of the perturbative and non-perturbative correc¬ 
tions. Especially, we will focus on how these exact expansions are obtained from various 
beautiful relations to Fermi gas system, topological string theory, integrable model and 
supergroup. 
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1 Introduction 


Chern-Simons theories play a central role in modern string theory. Two major subjects in 
early days are open string field theory (IJ and topological string theory |2j. It was known not 
only that Chern-Simons theory satisfies the axiom of topological field theory but also that 
the topological theory can be lifted to topological string theory. Also, the beautiful gauge 
invariant structure was used to construct the covariant open string field theory. Both of these 
developments are deeply related to understanding non-perturbative aspects of string theory. 

The relation between the Chern-Simons theory and the non-perturbative string theory 
becomes even more important in the supersymmetrization. It was known that up to Af = 3 
the supersymmetric Chern-Simons theory can be constructed for any gauge group and any 
representation M- Also, it was noted that when the levels are summed up to zero, the 
theory enjoys the conformal symmetry Earn. After some special arguments of the A/ = 4 
enhancements m, finally it was found that when the gauge group is U(N)kXU(N)_k with the 
subscripts denoting levels k and —k and the matters are in the bifundamental representation, 
the supersymmetry is enhanced to Af = 6 [10] , which is called the ABJM theory. Later, it was 
also found mm that even when the ranks of two factors in the gauge group are different, 
U(Ni)k x U(N 2 )_ k , the Chern-Simons theory still enjoys the enhanced supersymmetry Af = 6, 
which is called the ABJ theory. So we shall often refer to them inclusively as ABJ(M) theory. 

The discovery of this maximal supersymmetric conformal Chern-Simons theory enables us 
to open up another window to non-perturbative aspects of string theory. In the study of the 
non-perturbative effects in string theory, we were naturally led to a conjecture that there exists 
an eleven-dimensional theory, which is dubbed M-theory ra- It is expected that this theory 
is the mother theory with fundamental membrane excitation from which our perturbative 
string theory stems after compactifications, though only little is known about this mysterious 
theory. From the supergravity analysis, it is known that, besides the fundamental electric M2- 
branes, there exist solitonic magnetic M5-branes. Also, from the AdS/CFT correspondence, 
we know that the degrees of freedom of N M2-branes and N M5-branes are AC/ 2 and N 3 
respectively [14]. 

The ABJM theory describes a stack of N multiple M2-branes put on the target space 
C 4 /Zfc, while the ABJ theory describes a composite system of both Afi multiple M2-branes and 
N '2 — N± multiple fractional M2-branes for N 2 > N±. For other superconformal Chern-Simons 
theory with less supersymmetries [6HTHT5]. we expect that they describe multiple M2-branes 
on various backgrounds with less supersymmetries. It is then natural to ask whether we can 
reproduce the degrees of freedom by studying the free energy of the ABJM theory in the large 
N limit. Also, it is interesting to see how the large N corrections behave and whether there 
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is any interesting physical implication in it. 

In contrast to Chern-Simons theory, let us shortly comment on Yang-Mills theory. Along 
with the Chern-Simons theory, Yang-Mills theory also plays a crucial role in modern string 
theory. Especially, it is known that multiple D-branes are described by the maximally super- 
symmetric Yang-Mills theory with the gauge group U(N ) [TB]J- This was found by identifying 
the dynamical degrees of freedom of D-branes as open string excitations appearing in each 
pair of the multiple D-branes. The free energy N 2 in the large N limit is interpreted as the 
number of the matrix elements and matches with the supergravity computation [14] . 

Although the maximally supersymmetric Chern-Simons theory selects a special gauge 
group and a special representation, the maximally supersymmetric Yang-Mills theory is de¬ 
fined for general gauge groups and representations. Then, it is natural to expect that the 
gauge group and the matter contents of the ABJ(M) theory are special. In fact it was ob¬ 
served in HU that, when we look at the matrix model appearing after applying the localization 
techniques, the gauge group and the matter content have a good interpretation in the super¬ 
group U(Ni\N 2 ). Interestingly, the relation to the supergroup U(Ni\N 2 ) persists even in the 
one-point function of the half-BPS Wilson loop. The underlying relation to the supergroup 
may go beyond the accidental relation after the localization techniques and plays some roles 
in understanding the ABJ(M) theory itself. 

In the supersymmetric theory, it was found that when we compute the partition function 
or the correlation function preserving supersymmetries, the infinite-dimensional path integral 
reduces to a finite-dimensional matrix integration [183. For more careful explanation of the 
localization techniques, see a review [T9] in this volume. 

After the localization techniques, the partition function of the ABJ(M) theory is reduced to 
a matrix model m- Although many works are concentrated on the matrix model expression 
itself, solving the matrix model is an interesting and important subject. In physics solving one 
theory classically means finding the exact solutions to the equation of motion, while solving 
one theory on the quantum-mechanical level means performing the path integral exactly. Since 
the localization techniques already reduce the infinite-dimensional path integral into a matrix 
integration, solving the ABJ (M) theory amounts to performing the matrix integration exactly. 
This is our main results in [2TH24] and we shall explain it in more details in the remaining 
part of the introduction. 

Before explaining the details, let us try to discuss our motivation and expectation for the 
results. First of all, the full expression of the ABJ(M) partition function may provide some 
hidden structures of the M2-branes. After reproducing the predicted N 3 / 2 behavior in the 
large N limit from the matrix model [[25], it is natural to ask how the corrections look like. 
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Especially, compared with the matrix interpretation of the N 2 behavior for the D-branes, it 
is perplexing to find the iV 3 / 2 behavior for the M2-branes. We expect that finding out all the 
corrections will give us some hints to this mysterious N 3 / 2 behavior of the M2-branes. Prior 
to the discovery of the ABJ(M) theory, there were attempts in finding out the worldvolume 
theory of the M2-branes from a generalization of the Lie algebra [2B1I27] . Though it turned 
out the algebra is very restricted and is not explored extensively any more, these works may 
suggest a hidden novel structure in M2-branes to be discovered. 

Secondly, the ABJ(M) theory is one of the maximally supersymmetric theories. If we 
consider that the uniqueness of string theory stems from the maximum in the supersym- 
metrization, we expect that the maximally supersymmetric theories will play a crucial role 
in understanding string theory. At this point it should be important to study carefully all of 
fundamental quantities in the maximally supersymmetric theories. 

Thirdly, let us stress the importance of the exact solution. The main method in studying 
physics is the perturbation theory. However, in a standard situation, the physics beyond the 
perturbation theory has a very rich structure. We expect that the exact instanton expansion 
will provide us a window to understand non-perturbative effects. 

In the remaining part of the introduction, let us shortly recapitulate the main results of 
a series of our works [2T1 [24] . We will see that the results provide some insights to the above 
questions. 


1.1 ABJ(M) matrix model 


After applying the localization techniques, the partition function of the ABJ(M) theory reduces 
to a matrix model 


Z k (N 1 ,N 2 ) = 


- 1 ) 


ijVi(Vi — 1 )+±N 2 ( 7 V 2 — 1 ) 


Ir n i + jv 2 i= d 


Ni\N 2 \ 

JVi n 2 

n d ^ n d, a 

3 =1 
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( 1 . 1 ) 


where n& and n£ denotes the product taking over all of pairs (i,j) satisfying 1 < % < 
j < Ni and 1 < i < j < N 2 respectively, while the product nfj’^ 2 denotes n^i 11^= i • The 
integration is defined by 
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The derivation of this partition function is originally done in [20] • Here we shall not go into 
the details of the derivation and instead take this ABJ(M) matrix model (II.ip as our starting 
point. 


1.2 Supergroup structure 


This matrix model comes from the maximally superconformal Chern-Simons theory. Let us 
shortly comment on the beauty of this matrix model H7j. In (II.ip the most complicated part 
of this partition function is probably the measure 

n y,.(2 sinh f n%(2 sinh ^ f 

n; v J i,JV2 (2cosh^zki) 2 ' 1 ' ’ 

This measure can be simplified largely by dropping one factor of the gauge group 

At _ 2 

n(2 S inh^) . (1.4) 

i<j 

This measure is the famous one of the Chern-Simons matrix model. If the reader is not familiar 
with it, she can replace the hyperbolic function by a rational function 

At 

n<*-*)’• (i-5) 

i<j 


This is nothing but the U(N { ) invariant Vandermonde measure appearing in the diagonal- 
ization of the At x At Hermite matrix in the Gaussian matrix model. We can restore the 
hyperbolic function by exponentiating the variables /i into e ±At . If we want to restore the de¬ 
pendence of the variables v, all we have to do is the supersymmetrization. Namely we replace 
the invariant measure of U(N\ ) by that of the supergroup U(Ni\N 2 ) 

n ,"> 4 - ft ) 2 n * 


n<j 


Vi 


nTAft 


(1.6) 


L h3 


Vi 


After exponentiating /i, : , is 3 into e ±Mi , e ±Uj respectively, we come back to the original measure 
of the ABJ(M) matrix model. The explanation is summarized in table CD 


In this sense, we can say that the ABJ(M) matrix model is the generalization of the 
most fundamental Gaussian matrix model with the simultaneous deformation of the super¬ 
symmetrization and the Chern-Simons deformation. As a message, we would like to stress 
that, aside from the physical interpretation, even purely from a mathematical viewpoint, the 
ABJ(M) matrix model is one of the most fundamental matrix models which deserve careful 
study. 
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Table 1: The measure of the ABJ(M) matrix model can be regarded as a simultaneous de¬ 
formation of the supersymmetry deformation and the Chern-Simons deformation from the 
standard Vandermonde measure. 

1.3 Complete large N expansion 

Now let us summarize the result found in a series of works [2TH251 [28H39] for both the 
ABJM case and the ABJ case. Let us reparametrize the arguments of the partition function 
Zfc(Ai, N 2 ) in (11.11) as N) with (N\, N 2 ) = (N,N + M) and define the grand canonical 
partition function with 

OO 

E kiM { f i) = ^2e> tN \Z k7M {N)\, (1.7) 

N =0 

and the (modified) grand potential J/ C) m(A 4 ) with (22J 

OO 

ZkMri = X) e J ^ + ^ in \ (1.8) 

n =—oo 

by regarding the rank N as the number of particles and introducing the chemical potential /i 
dual to N. Note that the phase dependence of the partition function Zk,M(N) was studied 
carefully in |20j and to make contact with the topological string theory we have to take the 
absolute value of the partition function in defining the grand potential [2Hj. It is important 
to notice that the large N physics in the canonical ensemble is governed by the large fi limit 
in the grand canonical ensemble, thus below we focus on this expansion. 

Note that the grand canonical partition function defined originally in (II. 7ft is symmetric 
under the 27 ri shift of the chemical potential 

/j, —> n + 2ni. (1.9) 
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Due to the periodicity the grand partition function contains an oscillatory behavior in an 
apparent sense [22]. The (modified) grand potential in (ll.8[) is defined by getting rid of this 
oscillation. In the inverse transformation, the partition function is given in terms of the grand 
canonical partition function ’Ek,M{fi>) by 

\Zk,M(N)\ = / —^ 7 Sfc ! M(A i )e _MAr , (1.10) 

J-iri 

while in the grand potential (1 1.12 j) by 

\Z k , M {N)\ = (1.11) 

J-ooi 2 m 

This difference plays crucial roles, for example, in discussing the generalizations of the ABJM 
theory in [00]. 

Let us first display our main result, and explain how it was found and what the physical 
implication is later in the subsequent subsections. The large /i results of the grand potential 
can be summarized as follows. The grand potential is split into the “perturbative” part and 
the “non-perturbative” part 


vw = •4 p ,«v s )+•yw 1 )- 


( 1 . 12 ) 


where “effective” chemical potential /i cS will be explained later. The perturbative part is given 
by 


^k,M (h) — _ 3 _ ^ 3 ^ k ’ M ^ + 

where the coefficients Ck and B^m are 
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For integral k, the integration is performed explicitly. 
The non-perturbative part is given by 

1 d 


r’’(v) = F ,op (T\g,) + 


2m dg s 


s,c NS (-;- 

9s 9s 
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Here the functions F top and F NS are defined by the limits 


F top (T; r) = lim F ref (T; n , r 2 ), 

T\^fT,T2^~T 

F ns (T; r) = lim 2^r 2 F ref (T; n, r 2 ), 

ri—>- t,7*2—>■() 
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F“ f (T;T,,T 2 )=£££jV* 
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q t =e Mn , q 2 = e 2 "^, q L =e l{n ~ n \ q R = e"' T ‘ +T “>. 


( 1 . 20 ) 


The functions in (11.171) have their origin in topological string theory. Namely, the function 
F rcf (T; Ti, r 2 ) is the free energy of the refined topological string theory, while F top (T ; r) and 
F ns (T ; r) are the free energy of the original topological string theory and the free energy in 
the so-called Nekrasov-Shatashvili (NS) limit. The integers are the BPS indices, which 

are topological invariants counting modes on the manifolds. The relation between two Kahler 
parameters T and the chemical potential /i is given by 




,1 M\ 


( 1 . 21 ) 


while the string coupling g s is identified to the level k as 


g s = 2/k. 


( 1 . 22 ) 


Note that, in (I1.12p . we need to substitute for the effective chemical potential 

/i — 2(—l)^~ M e~ 2/ NF 3 (l, 1, |, |; 2, 2, 2; 16(—l)t~ M e _2/i ), for even k, 


= 


p + e-%F 3 (l, 1, §, §; 2, 2, 2;-16e- 4 ^), 


for odd k. 


(1.23) 


It is interesting to find that the ABJ(M) matrix model has a deep relation to the topological 
string theory 03. 

This result, with the effective chemical potential plugged in, shows that the non-perturbative 
effects consist of terms 

OO 

r p (p) = (1.24) 

m,i =0 

0 , 0 ) 
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with fm.til- 1 ) being polynomials of /i. Here the effect of e~T is identified with the worldsheet 
instanton, which is a fundamental string wrapping the holomorphic cycle CP 1 in C 4 /Zfc on 
the gravity side [251 IT3] . and that of e~ 2fl is the membrane instantons, which is a D2-brane 
wrapping the Lagrangian submanifold MP 3 [29] • In particular, the coefficients ,/ m y(/i) with 
£ = 0 or m = 0 take the following form [52] 

fm, o(a0 = d m (k), foM = a e (k)fi 2 + b^k)^ + c £ (k). (1.25) 

Besides the worldsheet instantons and the membrane instantons, the full instanton series (11.241) 
contains the bound states of these two instantons. 

1.4 History 

Here we shall explain how the result (1 1.12 [1 with (jl, 13j) and (11.161) was found in chronological 
order. From the AdS/CFT correspondence, the N 3 ^ 2 behavior of the degrees of freedom of 
the N multiple M2-branes [T4j in the large N limit was known from the gravity side about 
two decades ago. The study from the gauge theory side to reproduce the N 3 / 2 behavior 
and investigate the corrections, on the other hand, was started quite recently [25] after the 
discovery of the worldvolume theory of the multiple M2-branes pjQ]. For the explanation, we 
shall mainly concentrate on the ABJM case. 

Since the partition function was reduced to a matrix integration from the localizations of 
the supersymmetric theories [2011441145] . we can apply’t Hooft expansion, the standard matrix 
model technique to it. The ’t Hooft expansion is to consider the ’t Hooft limit 

N —* oo, A = N/k : fixed, (1.26) 

and keep track of the 1/N corrections. Then, the free energy F(N,X) = logZfc j o(A r ) is given 
by 


F(N, A) = N 2 F 0 ( A) + Fi(A) + N~ 2 F 2 { A) + • • • , (1.27) 

which can be regarded as the genus expansion. From this viewpoint, it is not difficult to 
imagine that the’t Hooft expansion has a good interpretation in string theory. 

The first important work of the ABJM matrix model is the reproduction of the N 3 / 2 
behavior [25] in the ’t Hooft limit. In [251129j it was further noted that the expansion has 
infinitely many 1/N corrections 

F 0 (\) = \- 1 *+O(e -^), 
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(1.28) 


Fi(A) = Ai + logA + 0(e-^), 
F 2 (\) = \l + \ 2 + + Oie-^), 


where we only present the result schematically. Especially, the coefficients are not displayed 
and the ’t Hooft coupling constant is slightly shifted from the original one, but we will only 
come back to these points later. The study of the 1/N corrections is most easily done by 
observing the relation between the ABJ(M) matrix model and the topological string theory 
on local P 1 x P 1 [42] and studying with the holomorphic anomaly equation [46] from the 
topological string theory. Here it is interesting to notice a triangular structure in the expansion 
(ll.28jh Namely, by going to one higher genus from F g ( A) to F s+ i(A) we always gain only one 
more perturbative term and the perturbative terms in F g+ i{\) / F g {\) range between A 2 and 
A^ in the large A limit. This means that the large N expansion can be alternatively obtained 
by taking the large N limit with X 2 N~ 2a (0 < a < 1) fixedE If we look at the coefficients 
carefully, it is not difficult to observe that after introducing m 


A ren = A- 

24 


A 2 


3iY 2 ’ 


(1.29) 


the infinite expansion can be taken care of partially. So the remaining question is whether 
we can perform the other summation by dropping the non-perturbative terms and solving 
the perturbative holomorphic anomaly equation. In [30] it was found that if we drop the 
non-perturbative corrections all of the perturbative expansions can be summed up to the Airy 
function 


zr‘w = e A "C-‘ Ai <V(JV - B „, 0 ) 


(1.30) 


In [30j the constant Ak was not taken care of but this constant was later pointed out in [33] 
and studied carefully in [TTj. There are many arguments on the meaning of the Airy function 
appearing in the summation of the perturbative correct ions 0 Technically it suggests the 
viewpoint of the grand canonical ensemble. Namely, the integral representation of the Airy 
function is 

Ai(N) = / — e^ 3 -^, (1.31) 

Jc 2 ni 

where the integral contour C comes from the infinity in the region with angle —7 t/2 < 6 < 
—7t/ 3 seen from the origin and goes to the infinity in the region with angle 7 t/3 < 6 < n/2. 
Comparing this expression with the inverse transformation (ll.lljh it is not difficult to find 

*See [17] for more discussions on the limits. 

^See a recent work [48] for interesting generalizations and discussions. 
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that the complicated Airy function is mapped to a simple cubic function (11.13j) . This suggests 
that the full expression is clearer in the grand canonical ensemble. 

As pointed out in [28], there are dissatisfactions with the ’t Hooft expansion. In the ’t 
Hooft expansion, the large N limit is taken by fixing the ’t Hooft coupling A = N/k. However, 
in approaching to the M-theory regime in the large N limit, the parameter we want to fix is 
the level k itself which characterizes the M-theory background C 4 /Z*., 

N —> oo, k : fixed. (1-32) 


These ideas were incarnated in [32], by rewriting the partition function of the ABJM theory 
with gauge group U(N) k x U(N)_ k into that of a Fermi gas system with N non-interacting 
particles whose Hamiltonian is given by 


111 
y/2 cosh § 2 cosh | cosh § 


(1.33) 


with the canonical commutation relation [q,p] = ih where the Chern-Simons level plays the 
role of the Planck constant, h = 2nk. 


This beautiful rewriting directly suggests another analysis for the study of the ABJM 
matrix model. Namely, we can study the grand canonical ensemble of the ABJM matrix 
model from the WKB expansion around h = 0 1 32] . In this expansion, we first expand the 
whole grand potential in the small k region 

JM = + kJ {2 \p) + fc 3 J (4) (/i) + • • • , (1.34) 

and further study each term separately in the large /i expansion. After that, we collect all 
terms of the same order in /i and hope that we can sum up the k expansion. In this sense we 
can separate the small k expansion and the large /i expansion and take the large /i limit by 
keeping the level k fixed. 


Both the ’t Hooft expansion and the WKB expansion are powerful to study the matrix 
model to some extent. On one hand, the ’t Hooft expansion detects not only the perturbative 
part of the Airy function but also the worldsheet instanton part, though is sensible to neither 
the membrane instanton part nor the bound state part. This is because after plugging in the 
stationary condition for /i 

N = C k fi 2 , (1.35) 

0)1 

the worldsheet instanton factor and the membrane instanton factor in (jl.24j) become 

e-t-e" 2 ^, (1.36) 
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On the other hand, the WKB expansion detect only the perturbative part and the membrane 
instanton part, but not the worldsheet instanton part or the bound state part. So, unfortu¬ 
nately, for the bound state of the worldsheet instantons and the membrane instantons, neither 
the ’t Hooft expansion nor the WKB expansion is applicable. The only method we have is 
the numerical analysis and the pole cancellation mechanism, which we will explain in the next 
subsection. 


1.5 Pole cancellation mechanism 


In determining the non-perturbative effects, the so-called pole cancellation mechanism plays 
an important role. Let us recapitulate it shortly. For this purpose we display the explicit 
forms of the first few worldsheet instantons in (II.25ft 


di(k) = 


Sill 


2 2tt ’ 


d 2 {k) = 


2 sin 


2 4tt 


Sill 


2 2t r ’ 


and membrane instantons 


ai (k) 
a 2 {k) 


4 cos i y 
n 2 k ’ 

8 + 10 cos n k 
n 2 k 


bi(k) 
b 2 {k ) 


2 cos 2 f 
7T sill ' 


4(1 + cos 7T k) ^17 + 24 cos tt k 


7T 2 k 


2n sin7r k 


ci(k) 

c 2 (k) 


(1.37) 


J 


(1.38) 


It is not difficult to fold that the coefficients of the m-th worldsheet instanton contain poles 
when e Z, while the coefficients of the f'-th membrane instanton contain poles when 
2 £k e Z. However, it is interesting to observe that whenever one coefficient is divergent the 
other coefficient is always divergent as well and totally the coefficient is completely finite. This 
property was first observed in [22] and the importance was stressed in [35], 

This property is probably due to the supersymmetries. Since there is no phase transition 
in theories with enough supersymmetries, the partition function should be finite. From the 
viewpoint of applications, this property can be used to fix first few instanton coefficients. In 
fact the first few coefficients of the membrane instantons listed above were found from this pole 
cancellation mechanism. The existence of the membrane instantons was first observed in [49]. 
It is surprising that the first quantitative study of the membrane instanton is performed in 
this way. 

We can proceed to the bound states and find a similar pole cancellation mechanism there. 
Namely, instead of the simple pole cancellation mechanism between the worldsheet instanton 
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and the membrane instanton, for higher instanton effects the bound state of the worldsheet 
instanton and the membrane instanton appears as well. Hence, we need to consider the 
cancellation among all of the bound states besides the pure worldsheet instanton and the pure 
membrane instanton. 

In |23] we find that the effect of the bound states can be incorporated into that of the 
worldsheet instantons by shifting the chemical potential as in (jl.23j) . As a bonus, after shifting 
the chemical potential n into the effective one /j, e g, we find that the quadratic polynomial 
coefficients of the pure membrane instanton fo/(n) fl 1.2511 become linear polynomials, where 
the coefficients of the linear terms and constant terms obey a derivative relation. Finally, 
in EH, all of these structure are reexpressed in terms of the refined topological string theory 
as in (11.121) . 

I. 6 Contents 

The main purpose of this paper is to review the beautiful aspects of the ABJ(M) matrix model 
and its generalizations. After reviewing the history of the progress in section 11.41 and section 

II. 51 it is not difficult to imagine that the explanation would be the simplest if we start with 
reviewing the Fermi gas formalism. We shall do this in the next section first for the ABJM 
matrix model and then generalize to the ABJ matrix model. After that, in section [3l we shall 
proceed to various analysis including the WKB expansion, the exact evaluation of the partition 
function, the numerical study of the grand potential, the pole cancellation mechanism and the 
viewpoint from the topological string theory. Finally, in section [4] we summarize some recent 
progress in the supersymmetry enhanced case and more general cases. 


2 Fermi gas formalism 

In this section, let us explain how the Fermi gas formalism is applied to the ABJM and ABJ 
matrix model. 

2.1 ABJM matrix model 

In this subsection, we shall explain the Fermi gas formalism for the ABJM matrix model [32j. 
Namely, we shall study (11.111 by restricting ourselves to the case of = N\ = N and define 

Z k (N) = Z k (N,N). (2.1) 
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Using the Cauchy identity 


det 

Ui + Vj 

one can show that (11.11) is rewritten as 

D n h D n v 


IW«i - u i)( v i - v i) 


Z k {N ) = 


N\ N\ 


det 


YliAUi + Vj 


det 


( 2 . 2 ) 


2 cosh — - Uj J l<i<N 

2 1 <j<N 


2 cosh — / i<*<(v 

2 i<j<iv 


(2.3) 


Note that the two determinant factors are identical, though we write them separately, so 
that the structure is clearer. Using a determinantal formula in [28], we can combine the two 
determinants 


Z k (N) = 


D N p 
N\ 


det 


2 cosh 2 cosh hAh J i<i<N 

2 2 1 <j<lV 


(2.4) 


Note that here o stands for the integration 

<p o ip = / Dv(j){y)il){y). 


(2.5) 


In other words, we regard two ingredients (2 cosh ) 1 and (2 cosh Up) 1 as matrices with 
the continuous indices /i and za As in the discrete case, where we perform the matrix multi¬ 
plication by contracting the indices, here the contractions of the continuous indices are done 
by the integrations given in (11.21) . 

After expanding the determinant of (12.4[) and rescaling the integration variables p — x/k, 
v = y/k, we find that it can be expressed as a sum over permutations Sn 

N 


, -1 


Zk ^ N) = m S (_1) / 

ctESn ^ i =1 


(2.6) 


with p(x,x') given by 


, /\ /' dy _j_ x 2 1 1 

'<*■*) = hi e 


i „.2 


g 47r/e y — 


2fc 


2 cosh - 


2fc 


If we use the Fourier transformation 

1 1 


k 2 cosh 


= (£ 


2 cosh | 


i»>, 


(2.7) 


( 2 , 8 ) 


by introducing the coordinate operator q and the momentum operator p satisfying [q. p] = ih 
with h = 2nk and the coordinate eigenstate 


q\x) = x\x), 


(2.9) 
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normalized by {x\x') = 2 / n8{x — x'), we find that p(x,x') can be expressed as 

1 i *2 1 


p(x,x') = (x|e 2 " 9 


2 cosh 2 


— —a^ 
-e 2hi . 


2 cosh | 


| a/). 


After some further canonical transformations with the help of the formula 
e irf(p)e-i r = f(p-q), g(q)e~A’ l ‘ =g{q+p). 


( 2 . 10 ) 


( 2 . 11 ) 


and a similarity transformation, we finally find that p(x, y ) can be interpreted as a density 
matrix of a quantum mechanical system 

p(x,y) = (x\e~ H \y), (2.12) 


with (11.331) . which leads to the expression 


p(x,y) = t- 


(2.13) 


k 2 cosh | 2 cosh ^J2 cosh 2 

The formula (12.6p can be regarded as the (—l)°-graded partition function of a particle prop¬ 
agating from a state labeled by the coordinate Xi to the permutated coordinate x^. This is 
nothing but the partition function of a Fermi gas system. To avoid unnecessary complexities, 
hereafter we shall often drop the hats for the operators q and p. 


ft is convenient to define the grand partition function S*.(/i) by introducing the chemical 
potential p, conjugate to N and summing over N. It turns out that Sfc(p) can be written as 
a Fredholm determinant of p 

OO 

ZM = J2 eN " Z ^ N ) = det (! + e V), (2.14) 

N =0 

and the grand potential Jk(p) = logSfc(/i) is given byEl 

1=1 

The expression of the Fredholm determinant (12. 14ft and its expansion (12 .1 5H can be shown 
by relating the coefficients of J^^^Tr p £ ) me , to the combinatorics of distributing N into an 
assembly of mi cycles of length £ satisfying = ^ |5Q] . 

In terms of the eigenvalues E n (n = 0,1, • • •) of the Hamiltonian H, the grand partition 
function is given by 

OO 

= IF 1 + el, ~ E ")’ < 2 - 16 ) 

n= 0 

*Note that this grand potential Jk{y) is different from the modified one defined in (11.81) . 
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hence the problem boils down to the diagonalization of the Hamiltonian. The study of the 
spectrum of the ABJM Hamiltonian was initiated in |2T] using the Hankel matrix representa¬ 
tion of p, which was later generalized to other models. 

Note also that the partition function can also be rewritten as 



(2.17) 


which is useful to discuss the generalization to the ABJ model. 

2.2 ABJ matrix model 

In the previous subsection, we have given the Fermi gas formalism for the ABJM matrix model. 
Using this Fermi gas formalism, we shall explain how the exact instanton expansion was found 
in the next section. Before it, let us comment that the Fermi gas formalism is also applicable 
to the ABJ matrix model where the gauge group U(Ni)k x U(N 2 )-k has different ranks. In 
generalizing the discussions in the previous subsection, there are two Fermi gas formalisms 
for the ABJ matrix model. One is to dress the density matrix in (12.131) by extra factors with 
the traces kept unchanged, while the other is to keep the density matrix and introduce the 
endpoints for the density matrix. 

Since the extra physical ingredients of the ABJ theory compared with the ABJM theory 
is the fractional M2-bra.nes, two formalisms amount to expressing the fractional branes by 
changing the closed string backgrounds or from the open string excitations. Hence, it is 
reasonable to call these two formalisms as closed string formalism and open string formalism 
respectively. 

After the perturbative coefficients C*., (1 1.141) and A & (jl. 15j) and the worldsheet in- 

stantons were found from the consistency with the ’t Hooft expansion [[25], the topological 
string on local P 1 x P 1 and the ABJM limit M —> 0, we can use either of these Fermi gas 
formalisms to further study the membrane instantons and the bound states. 

2.2.1 Closed string formalism 

Let us start the study of the ABJ matrix model (11.11) with different ranks. In what follows, 
we will set Ni = N, N 2 = N + M and assume k, M > 0 without loss of generality. First let 
us note that the partition function for IV = 0 is nothing but the bosonic U(M) Chern-Simons 
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matrix model, which is given by |51j 


Z k {Q,M) 


M 

2 


M—l 

n 

S=1 


2 si ”t) 


Ms 


(2.18) 


As discussed in [MISZIES], if we define Z k (N, N + M) to be the absolute value of the partition 
function Z k (N,N + M) normalized by the Chern-Simons case (N = 0), 


Z k (N,N + M) 


Z k (N, N + M) 
Z k (0, M) 


(2.19) 


we find that Z k (N, N + M) is given by an integration 


N 


1 [' d N r -n- 

Z k (N, N+M) = — I J] tank 


N 


2 Xi Xj 


N\ J (4? rk) N 


i<j 


2k 




( 2 . 20 ) 


i= 1 


with a dressing factor 


M — 1 
2 


V{x) = 


ea 


-l) M e 


tanh 


x + 2n is 
2k ’ 


( 2 . 21 ) 


which is very similar to (12.1 7ft in the ABJM matrix model. The Seiberg-like duality of the 
ABJ theory 


U(N) k x U(N + M)_ k U(N)_ k x U(N + k-M) k (2.22) 

implies that Z k (N,N + M) is invariant under the exchange M •<->■ k — M with k and N 
fixed. Also, note that the original partition function Z k (N,N + M) (as well as the Chern- 
Simons case Z k ( 0, M)) vanishes when k < M, which is consistent with the conjecture that the 
supersymmetry of the ABJ theory is spontaneously broken for k < M. 

Again, using the Cauchy identity, Z k (N, N + M) can be recast as a Fermi gas form 

1 r rl N N 

Z t (N,N + M) = —J2(-ir / < 2 - 23 ) 

■ o-eSjv J i= 1 

with the density matrix (j2.13j) changed into 

P(x,y) = lV v (x) 2c ^ hE _y V v (v)- (2-24) 

Using this dressed density matrix, the grand partition function E kjM (/j) in (II.7ft is given by 

5 fe ,iu(^) = |^fc(0, M)\ det(l + e^p). (2.25) 
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2.2.2 Open string formalism 


Having seen the description of the fractional branes by changing the closed string background 
in the previous subsubsection, here let us turn to the description from the open string |38|. 
This formalism has a benefit in the compatibility with the vacuum expectation values of the 
half-BPS Wilson loop, since we end up with an ultimate Fermi gas formalism which includes 
both the M2 fractional branes and the half-BPS Wilson loop. As a corollary of this formalism, 
we can reproduce, on one hand, the Fermi gas formalism for the vacuum expectation values 
of the half-BPS Wilson loop [52], and on the other hand, that for the partition function of the 
M2 fractional branes [ 38] . 

To state the results we need some preparations. Instead of the partition function itself, 
since the formalism incorporates the insertion of the half-BPS Wilson loop as well, let us 
consider the vacuum expectation value of the half-BPS Wilson loop. After the localization of 
the supersymmetric theory, it was found [171 [20] that the insertion means to include a loop 
operator in the matrix model 


(s Y )k(N 1: N 2 ) = 


(-l^JW-lJ^JVa-r) r d N lfi d N 2p 


N X \N 2 \ 


X 


n 


J ( 271)^1 ( 27^2 

N \ 2 sinli ^ 2 sinh " n 2 


s y (e'V) 


i<] 




nr=im=i2cosh 


f-Lj Vj 
2 


1 ii 2 —v 2 v 2 \ 
£ 4 tt \z^i=l H'i 2^j=\ j ) ' 


(2.26) 


Here the loop operator sy(e /i |e I/ ) = Sy(e w , e M2 , • • • , e MJV i \e Ul , e 1 ' 2 , ■ ■ ■ , e VN *) is the so-called su¬ 
persymmetric Schur polynomial, which is nothing but the character of the hidden supergroup 
U (A51W). This is a symmetric polynomial with N\ variables /i,; and N 2 variables Vj, la¬ 
beled by a Young tableaux Y. It is probably the simplest to find out this supersymmetric 
Schur polynomial by expressing the original Schur polynomial in terms of a diagonal matrix 
diag(e w )^ 1 and the trace tr (in other words, the power sum symmetric polynomials) and re¬ 
placing the diagonal matrix by a diagonal supermatrix diag((e Mi )^ = 1 1 |(— e Uj )^ 1 ) and the trace 
by a supertrace. For example, let us consider the supersymmetric Schur polynomial of the 
second symmetric tensor in the supergroup U(Ni\N 2 ) with N { and N 2 being large enough 


s m (e M ) = ^(tr diag(e^)) 2 + ^ tr diag(e 2/i ). (2.27) 

Then after the replacement we find that the supersymmetric Schur polynomial can be decom¬ 
posed into the original one by 


s (e^l e v ) = s (e'H + s (e'Hs (e u ) + s (e u ). (2.28) 

m m □ □ p 
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Also, let us define the vacuum expectation value in the grand canonical ensemble 


MS, W = Y X + M)z N , (2.29) 

N =0 

as for the partition function. We further normalize it by the ABJM grand canonical partition 
function 


(sy) 


GC 

k,M 


[sr]Mf(*0 

MW ' 


(2.30) 


Then, we can prove 


MkM^) = det(H m (z))i<p< M +r, (2.31) 

1 <q<M+r 


where H pq (z) is given by 


H Pt q(z) 


E tp ( 1 + zQP ) l E_ M+q _ i, for 1 < q < M, 
zEi p (l + zQP)~ l QE aq _ M , fori < q - M < r. 


(2.32) 


We have to explain the meaning of each expression. First of all, we regard (Q)u,p and 
{P)p )V as matrices and (Ej) u as a vector with continuous indices /i and is. Again, as explained 
below (12.51) . we multiply matrices with continuous indices by the integrations (11.21) . So each 
expression is a scalar constructed by the inner product of two vectors with many matrices 
sandwiched in between. The explicit form of (Q)^, (P) MjI/ and ( Ej) v is 


(QW = 


2 cosh ^ ’ 


{P)w = 


2 cosh ’ 


(E,)„ = e^. 


(2.33) 


Now let us turn to the integers l p and a q . In short, these integers are p-th leg length and q-Wi 
arm length, which appear in the Frobenius symbol (aia 2 ■ ■ ■ a r ; lil 2 ■ ■ ■ I t +m)- The Frobenius 
symbol for the ABJM U(N\N) case is found by drawing the diagonal line for the Young 
diagram and counting the number of boxes horizontally or vertically for the arm length or the 
leg length, respectively. For the general ABJ case with U(N\N + M), all we have to do is to 
shift the diagonal line by M boxes. Note that, after shifting the diagonal line, the number of 
the arm lengths and that of the leg lengths are not equal any more. Hence, we cannot consider 
the determinant of the square matrix. We can restore, however, the equality of the numbers 
by consider the “inverse” arm length as well. The subscript — M + q — 1 appearing in (12.32ft is 
nothing but this inverse arm length. Therefore, we can say that the lower expression in (12.32(1 
comes from the usual arm length while the upper one is from the inverse arm length. (See 
e.g. [38,52] for figures explaining the Frobenuis symbol.) 
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Before explaining how the expression is obtained, let us first see how this formula is beauti¬ 
ful and useful simultaneously. First let us consider the ABJM case N 2 = N\ with a non-trivial 
insertion of the loop operator. Then, we find that M — 0 and the upper expression of H p q is 
absent. If we further use the same formula for the hook Young tableaux 

( s K;/ p ))m = zE lp (l + zQP)~ l QE aq , (2.34) 

then, as a corollary of the above formula, we can express the vacuum expectation value of 
a general loop operator as the determinant of those of the hook loop operators constructed 
by taking all combinations of the arm lengths and the leg lengths from the original Young 
tableaux [52] 

MS ( 2 ) = det (( s KY))£o(-))i<p<r- (2-35) 

l<g<r 

Note that, without taking the vacuum expectation value, it is a classical mathematical result 
stating that the character itself satisfies the same type of formula 

S Y = det (s(a,;Z p )) l<p<r, (2.36) 

l<q<r 

which is called Giambelli formula. Hence, in other words, what we have shown is that we 
can put to each loop operator in the Giambelli formula (I2.36P the normalized vacuum expec¬ 
tation values in the sense of the grand canonical ensemble. Namely, the normalized vacuum 
expectation values defined in (12.301) is Giambelli compatible. 

Next, let us turn to another special case of the partition function with no insertions sy = 1. 
Then, we find 

(!)fc M) = det(£ , M -p(l + zQP)~ 1 E_ M+q _ l ) l < p <Mi (2.37) 

i <g<M 

since l p = M—p in this case. If we regard the trace operator without endpoints as closed strings 
and the meson operators EQPQP ■ ■ ■ E which appears by expanding the expression in (12.32ft 
as open strings, this formula seems to state that the fractional M2-branes can be constructed 
from the open strings. This is reminiscent of the fact that classical D-brane solutions can be 
studied from open string field theory [53] . 

So far we have stressed the beauty of this formula. It is also interesting to note that this 
formula is actually very useful. In section [3721 we shall see that the special form of the density 
matrix is very useful to compute the powers of the density matrix. Our formula in (12.32ft also 
shows the same computability. Namely, due to the expression of the meson operator, we can 
start with a vector E and multiply matrices Q and P to it one after another without difficulty. 
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1 


1 

1 

1 \ / 

1 

1 

Ul+'Ul 

Ul+V 2 

U\+Vm 2 



Vl+Ul 

V2+U1 


Vn 2 + u 1 


\ 


det 


1 

UN-, +^1 

N 2 —N 1 —I 

U 1 

N2 — N1 — 2 
U 1 


1 

UN 1 +V 2 

„N 2 -N!-1 

u 2 

..N 2 -N 1-2 

u 2 


1 

un 1 +vn 2 
„N 2 -N!-1 

v n 2 

N 2 -Ni-2 

u n 2 


det 


1 

„,N 2 —N\—1 

U 1 

„.N 2 —Ni—2 

U 1 


1 

1>2+«JV 1 

N 2 —Ni — 1 

_y2_ 

Ni-N^-2 

u 2 


1 

vn 2 +u Nl 
n ,N 2 -N i-l 
V N 2 

n ,N 2 —N 1 —2 

U N 2 


Vi v 2 ■■■ Vn 2 


Vi v 2 ■■■ Vn 2 


\ 1 1 


1 / V 1 1 


1 


Tabic 2: An explanation how the EQPE term in the expansion of £'(1 + zQP ) 1 E in (12,32ft 
appears. 


Let us now turn to the explanation how this expression is obtained. First let us explain 
the formula for the partition function (1)*? m(^)- Again our starting point is the matrix model 
m- Instead of the Cauchy determinant formula (12.2)1 . let us consider the following formula 


/ 


det 


Hi 


\ 


1 / l<i<iVi 
l<j<N 2 

(„,N 2 —Ni— i\ 
l V i flKiKNi-Nx 

l<j<N 2 




n 2 


li<j 


i<r 


nrr 2 K+^-) 


(2.38) 


This is a natural combination of the Vandermonde determinant and the Cauchy determinant, 
which reduces to the Vandermonde determinant for the special case of N\ = 0, while reduces 
to the Cauchy determinant (12.2ft for N 2 = N\. This formula can be proved by starting with 
the Cauchy determinant (12.2[1 with matrix size N 2 x N 2 and sending v ^ 2 , ujv 2 -i, ■ ■ ■, UjVi+i to 
infinity one after another. 


Then, after plugging in u t = e Xi and Vj = e % , we can reproduce the measure of the 
ABJ matrix model. Let us see how we can put the measure into a Fredholm determinantal 
expression as in the case of the ABJM matrix model. Since the substitution can be done at 
any stage, for simplicity, let us stick to the variables Ui and Vj. 

Since in the partition function (II.ip we have two factors of this measure, we have two 
determinants. (See table[2j) So there are a lot of terms in expanding the determinant. Suppose 
we pick up one element from the lower ( N 2 — A r 1 ) x N 2 sub matrix of the first determinant, say 
v ^ 2 -n 1 - 1 . This means that, in this term, we are not allowed to pick up the other elements in 
the same row or in the same column. This element is associated to the integration 

of the variable v\. To finish the integration of we need to decide which element to pick up 
from the first column in the second determinant. 
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If we pick up one element from the lower (jV 2 — Ni) x N 2 submatrix of the second deter¬ 
minant, say v^ 2 ~ Nl ~ 2 , again this means that we do not pick up the others in the same row 
or column. Then, we do not have other elements depending on the variable V\ and we can 
complete the v\ integration. 

Instead, suppose in the second determinant we pick up another element from the upper 
Ni x N 2 submatrix, say (ui + Ui) _1 . Then, again there are no other elements associated with 
the variable v\ and we can integrate v\. But now we are associated with the integration of 
another variable u\. To finish the u± integration, we need to pick up one element from the first 
row of the first determinant, say ( U\ + v 2 )~ l . (Note that we can no more pick up the element 
(Ui + Ui) _1 from the first determinant, because we already pick up u^ 2- ^ 1-1 .) Then, we are 
associated with the integration of the other variable v 2 . 

Now, in the second determinant, if we choose from the lower (N 2 — Ah) x N 2 submatrix, 
then we finish the v 2 integration. But if we choose from the upper A r 1 x N 2 submatrix, we 
need to finish another integration of the variable u. 

We can continue this process and find a term consisting of a series of matrix multiplications, 
with two monomials of v on the two ends connected by products of ( v + u)~ l (u + u) -1 . 
After plugging in Ui = e Xi and Vj = e Vj , essentially v on the two ends turns into E while 
(v + u)~ l {u + v)~ l becomes QP. This can be regarded as an open string with Chan-Paton 
factors on the two ends. Of course, there are also contributions from the pure traces without 
ends. Now the only problem is to count the combinatorial factor correctly. The answer can 
be found more efficiently by preparing a mathematical formula for the matrix multiplications. 
This was done in [35] . 

For the vacuum expectation values of the loop operator, we can use another beautiful 
formula [54], stating that the supersymmetric Schur polynomial can be expressed as a ratio of 
two determinants, where the determinant in the denominator is nothing but the Vandermonde- 
Cauchy determinant we encounter in (j2.38j) . Hence, by combining the supersymmetric Schur 
polynomial with the measure, we find that one of the two determinants appearing in the 
measure of the partition function is replaced by another determinant. Hence, we can repeat 
the process explained in table [2] and find out the formula (12.311) . 


3 Non-perturbative corrections at large N 

After presenting the Fermi gas formalism for the AB J (M) matrix model in the previous section, 
we can now start our study of the large /i expansion for the grand potential. For simplicity, we 
will mainly focus on the ABJM case. Since the expression as the partition function of a Fermi 
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gas system already suggests a statistical mechanical method, we first start with the WKB h 
expansion. After that we turn to the study of the exact values of the partition function at 
finite N and see how these exact values lead to the numerical study of the grand partition 
function. Finally, we compare with the free energy of the topological string theory and see 
how the exact expression is found. 


3.1 WKB expansion 


As emphasized in [32J|55], the Fermi gas formalism enables us to study the M-theory regime 
(11.3211 . One useful way to study this regime is to perform the semi-classical (or WKB) expan¬ 
sion (ll.34p of the grand potential J fc (/i) = logS fc (/i) around k — 0 in the limit /i —> oo. 

For example, the leading order term in the WKB expansion is easily found by 

replacing Tr p e in (I2.15P by the classical phase space integral [32] 


J {0) W = £ 

i=i 


(- I )*" 1 

t 




dqdp f 1 1 \ 

2n \ 2 cosh | 2 cosh f ) 


/1 1 1 3 e 2/i \ 

T 3 2 V2’2’2’ ’2’l6 ) 


6 2 m 

87T 2 


4-C3 


1,1,1,1; 1,1,2; 

’ ’ ’ ’ 2 ’ 2 ’ ’ 16 J 


(3.1) 


From the large /i behavior of the hypergeometric functions in the above expression, 

4 (0) M = S + f + ^ 

+ -^(- 6 /i 2 + 6 /i + 6 - n 2 )e~ 2f * + ^(-36/i 2 - 66/1 + 25 - 6 vr 2 )e” ^ + 0(e~ 6/x ), (3.2) 

67T Z ZlT z 


one can see that the perturbative part of J^(/i) is indeed a cubic polynomial in / 1 . The higher 
order correction of the WKB expansion can be systematically computed by the method of 
Wigner transformation [32], or the semi-classical analysis of the thermodynamic Bethe ansatz 
(TBA) equation [35]. It turns out that the coefficient Ck and Bkp in (11.131) do not receive 
corrections higher than Oik 1 ) and 0(k 3 ), respectively. On the other hand, the constant A*, 
in (I1.13P receives corrections of all order in k. 

We can continue this process to higher and higher h corrections and reproduce the Taylor 
expansion of the functions a^(k) and bt(k) in (11.381) up to a certain large order of k. The 
coefficient of the membrane instanton was first obtained in this way [321135] , 
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Zi(l) = 

Zi(5) 
Z 2 (l) = 
M i) = 

Z 4 (l) = 

Z«( 1) = 


1 

4’ 


Zi(2) = 


167t’ 
26 + 20ir - 9ir 2 
40967T 2 ’ 


Z,(3) 

Zi(6) 


^ *(*> = 
78 - 121t r 2 + 36 t r 3 


10 — 7r 2 

1024tt 2 ’ 


147456tt 3 


1 

8’ 

1 

12 ’ 

1 

16’ 

1 

24’ 


^2(2) = 
3»(2) = 


32t r 2 ’ 

—3 + 7T 

487T 


^2(3) = 


10 — 7r 


^2(4) = 


24 - 32tt 2 + 3t r 4 


Z 3 (3) = 


512ir 2 ’ ' ' 49152ir 4 

9 + 1087T — 64\/3vr 


!V1 


^4(2) 

^e(2) 


-8 + tt 2 -8 - 32tt + 11tt 2 

512tt 2 ’ 4U “ 8192tt 2 

54 - 5t r 2 v _ 189 + 192^377 - 125 t r 2 

5184t r 2 ’ ’ ~ 186624tt 2 


Table 3: The first few exact values of the partition function of the ABJM matrix model. 
See [22] for more exact values. 


3.2 Exact values of partition function at finite N 


To fully study the instanton expansion, let us start the computations of the exact values of 
the partition function at finite N. The density matrix is of the form 


p(Qi,Q2) 


E(qi)E(q 2 ) 
M(qi) + M(q 2 )' 


(3.3) 


This type of the integration kernel is quite ubiquitous and appears from time to time in 
modern physics. In the context of light-cone string held theory [56], the so-called Neumann 
coefficients, representing the overlap between one string state and two string state, takes this 
form. From the integrability viewpoint, this form has a close relation to the TBA equation 
and was studied in E71 


Among others, this expression of the density matrix implies that we can compute the 
powers of the density matrix without much difficulty. To explain the computation, we first 
rewrite (13. 3 j) schematically as 


{M, p} = E <g) E, 


(3.4) 


if we regard p, M and E as a symmetric matrix, a diagonal matrix and a vector with con¬ 
tinuous indices respectively. Then, by studying the commutator [M, p n ] for even n or the 
anti-commutator {M,p n } for odd n, we fold 


n— 1 

P n (<?l,<?2) = J] 

m— 0 


C p m E)(q 1 )(p n - 1 - m E)(q 2 ) 
M(qi) - (-1 ) n M(q 2 ) 


(3.5) 
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J7 = 


J7 = 


4/i 2 + p, + 1/4 


7T Z 


e- 4/i + 


4/i^ + 2/i + 1 


7T Z 


e~ 2fl + 


52/i 2 + p/2 + 9/16 

Y 2 

52/i 2 + /i +9/4 


+ 2 


e- 8/i + C>(e- 12 ^), 


2vr 2 


+ 2 


3 -4 M 


C>(e 


-6/0 


Jo ip = -e — 2e 3^ + 
8 3 


J7 = e" p + 


J 6 np = i e -§/* - 2e-t'* + 

3 


4/i 2 + /i + 1/4 20 

Y 2 + ~9~ 


4/i 2 + 2/r + 1 

e -V + IV 3 ^ + 

52/i 2 + p + 9/4 o 

2vr 2 

3 

47r 2 


e - 4 M_^ e -^ + C )( e -^) ) 

9 


3 -4 M 


0(e 


-5 <u\ 


4/i 2 + 2/i + 1 20 

37t 2 9 


3 -2m _ ^ e -f^ + C>(e-f^). 
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Table 4: The first few non-perturbative terms of the grand potential of the ABJM matrix 
model obtained from the numerical study. See [22ll23lj for more terms. 


Note that, this formula shows that, instead of multipling the matrices directly, we can pick 
up a vector E and multiply p to it one after another, which is much easier. 

Using this method, in |22j we obtained the exact values of the partition function up to 
{k,N m ax) = (1,44), (2, 20), (3,18), (4,16), (6,14). See table [3] for a few examples. By now with 
some small technical progress, we believe that we can proceed to more exact values. See [39] 
for many exact values of the ABJ partition function. 


3.3 Grand potential 

After obtaining many exact values, we can now proceed to studying the exact coefficient of 
the total non-perturbative effects [22]. The results are given in table [4] Since our analysis 
contains some guesswork, we shall spell out the full details, so that the reader can judge the 
accuracy by herself. 

Let us pick up the case of k = 4 for concreteness. Suppose that we already know the 
expansion structure of the non-perturbative effects 

J77) = 7i e M + ( a 2 7 + + 72)e 2p + 736 3p + • • • • (3-6) 

By rewriting the large /1 expansion of the grand potential in terms of the large N expansion 
of the partition function, we fold 

Z 4 (N) = zf\N) + zf\N) + zf\N) + • • • , (3.7) 
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Z={1/16,1/512-1/(64*Pi~2),11/8192-1/(1024*Pi~2)-1/(256*Pi) ; 

(*Add as many exact values as possible*) 

$MaxExtraPrecision=300; 

Xi [w_]:=Exp[(gal) w+(al2 mu~2+be2 mu+ga2) w~2]; 

(*Add higher order terms w~3, w~4, ... similarly when necessary*) 
k = 4; c[k_]:= 2/(Pi~2 k); b[k_]:= 1/(3 k) + k/24; 
a[k_] := If[EvenQ[k], ~(Zeta[3]/(Pi~2 k)) \ 

- 2/k Sum[m (k/2 - m) Log[2 Sin[(2 Pi m)/k]], {m,1,k/2-l}], \ 

-(Zeta[3]/(8 Pi~2 k)) + k/4 Log[2] \ 

- 1/k Sum[(k + (-1)~m (2m- k))/4 (3k- (-l)~m (2m- k))/4 \ 

Log[2 Sin [(Pi m)/k]], {m,l,k-l}]]; 

JO = c[k]~(-l/3) AiryAi [c [k] ~ (-1/3) (n-b[k])]; 

J1 = CoefficientList[CoefficientList[Series[Xi[w],{w,0,2}],w][[2]],mu]. \ 
Table[(-l)~m D [c [k]~(-1/3) AiryAi[c[k]~(-1/3) (n + 4/k - b[k])],{n,m}], \ 
{m,0,0}] // Simplify; 

J2 = CoefficientList[CoefficientList[Series[Xi[w],{w,0,2}],w][[3]],mu]. \ 
Table[(-l)~m D [c [k]~(-1/3) AiryAi[c[k]~(-1/3) (n + 8/k - b[k])],{n,m}], \ 
{m,0,2}] // Simplify; 

(*Define higher instanton terms J3, J4, ... similarly when necessary*) 


Table 5: Mathematica code for studying the coefficients of non-perturbative effects. 


where each term is given by 


4 0) (iV) = e A4 C 4 1 Ai C 4 ^(N - B 4 ) 


zi 1] (N) = 7 ie A4 Q f Ai C 4 1 (N + 1 - B A ) 


Zf (N) = ( a 2 d 2 N - fod N + 72 + 7T ) # Ai C 4 # (N + 2 - B A ) 


7i 


(3.8) 


with B a = B A) o (11-141) . Note that the argument of the first non-perturbative effect Z 4 \N) is 
shifted by 1 because of the exponent factor e -M in J 4 p (n) (|3.6[) . In the second non-perturbative 
effect Z 4 \n), instead of the polynomial of /x, we replace /i by — c/v- Note also that the lower 
non-perturbative effects come in from the expansion of the exponential function in the grand 
canonical partition function. 


Now let us fit the exact values against this function form. We first prepare as many exact 
values as possible. In the following analysis we shall utilize all the values with A^ max = 16 
in |22j. As long as we know the command “FindFit” in Mathematica, it is by no means 
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difficult to obtain the fitted result. Still the reader may find it useful to directly apply the 
given program. So we also list the program in table [5] Note that, since we are only interested 
in integral k, we have used the expression for ag(k) in |41j . Then, if we ask Mathematica what 
is the best fit for 71 by 

FindFit[Log[N[Z,50]], a[k]+Log[Abs [J0+J1]],{gal},n] 

we find 


71 = 0.98721718103894311608972098476342715916438226064401. 

So we have 2% accuracy for 7 x. We are sure that the reader is not satisfied with it. For example, 
it is not possible to distinguish n from \/l0 with 2% accuracy. But it is not difficult, to roughly 
estimate the errors of this value. Although we have exact values for the partition function, 
since we neglect higher orders of the non-perturbative corrections in the grand potential, at 
each order we always have an error of e~ 4/x,/k . After plugging in the stationary condition 
N = dJ/dn = Ck/i 2 (11.351) . we find the errors 

e -4„A = e -2^ _ 10 -6 ) (3. 9 ) 

with k — 4 and N ~ 10. Since we have quadratic coefficients at the next order, the errors are 
actually much bigger than we expect. But anyway, we can improve our results by including 
higher instantons 

FindFit [Log[N[Z,50]] ,a[k]+Log[Abs [J0+J1+J2]] ,{gal,al2,be2,ga2}-,n] 

Now the result is much improved 

71 = 1.00000003159520570311264581986409771275067296902425. (3.10) 

We can continue this fitting to higher and higher orders. For example, for 

FindFit[Log[N[Z,50]],a[k]+Log[Abs[J0+J1+J2+J3+J4+J5] ], \ 

{gal,al2,be2,ga2,ga3,al4,be4,ga4,ga5},n] 

we find 


71 = 1.00000000000000000000000012119938655342513138232292. (3.11) 

However, if we go to the sixth non-perturbative corrections, the accuracy suddenly drops 

71 = 1.00000000000000000343882868029802074330254814069423. (3.12) 
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The reason is that although we only find roughly 25 digit accuracy, we have already used up 
the 50 digits in the exact values of partition function. Therefore, we improve by preparing 
100 digits from the beginning 

FindFit[Log[N[Z,100]],a[k]+Log[Abs[J0+J1+J2+J3+J4+J5+J6] ] , \ 

{gal,al2,be2,ga2,ga3,al4,be4,ga4,ga5,al6,be6,ga6},n] 

Then we find 

7 i = 0.999999999999999999999999999999999998548662756071218409269666•••. (3.13) 

When we extend to the eighth non-perturbative corrections, an error message appears again 

Failed to converge to the requested accuracy or precision within 100 iterations. 

But we can increase to 150 digits 

FindFit[Log[N[Z,150]],a[k]+Log[Abs[J0+J1+J2+J3+J4+J5+J6+J7+J8] ] , \ 

{gal,al2,be2,ga2,ga3,al4,be4,ga4,ga5,al6,be6,ga6,ga7,al8,be8,ga8},n] 

to avoid the error message and find 

7l = 0.999999999999999999999999999999999999999999999999999994165937•••. (3.14) 

Note that we should end our fitting at this point. We have only 16 data and if we fit up to 
the eighth non-perturbative corrections we use up all the data and it does not make sense to 
fit with more unknowns. We believe that we can go a little further to find more data with 
good machines. However, at present let’s be satisfied with more than 50 digit accuracy. After 
we are satisfied with the first coefficient, we can simply plug in the correct value and proceed 
to finding out coefficients of higher orders. 

Here we pick up the k = 4 case and explain our numerical study and the guesswork in 
details. We hope we have convinced the reader with the high accuracy we can reach at this 
point. 


3.4 Comparison with topological strings 

Note that the exponents of the grand potential obtained in table 0] have exactly the same one 
as the worldsheet instanton e ~in ('ll .24j) . Hence we can follow the proposal from [42] to 




compare the result with the free energy of the topological string theory on local P 1 x P 1 which 
describes the worldsheet instantons. We find that the numerical values of the first instanton 
for k = 3,4,6 and the values of the second instanton for k = 6 match with the prediction 
in (11.3711 coming from F top in (11.161) . However, if we trust this expression and extrapolate 
to smaller k , we find that the match does not hold any more. Besides, the expression is 
divergent at k — 1, 2 for the first instanton and k = 1, 2,4 for the second instanton. In fact 
this divergence gives us an important clue to understand the whole instanton expansion. The 
exponents of the terms containing divergences are e _2(?Al , which are nothing but the exponents 
of the membrane instantons. Besides, if we look at the behavior of the worldsheet instantons 
around the divergence, we find 


limdi(k)e k 

fc — >2 


lim d 2 (k)e k 

fc—>4 


4 4(p + l) 2p 2 + 2/i +1 11 

7r 2 (/c — 2) 2 7 r 2 (k — 2) 7T 2 3 

8 _ 4(/i + 1) _ 2/i 2 + 2p + 1 _ 7 

7T 2 (/c — 4) 2 7T 2 (/c — 4) 27T 2 6 


(3.15) 


whose finite part looks already close to the expression of the numerical values. In this way, we 
can expect that the membrane is also divergent at these points so that totally the divergences 
cancel among themselves. From table |4] we also observe that the divergence of the membrane 
instanton e~ 2£/i for an odd integer i comes from even integers fc, while the divergence for an 
even integer i comes from all integers k. Also, the coefficients of the membrane instanton 
e - 2 £/i f or oc pj £ vanish at odd integers k. Analytically continuing to non-integral fc, from this 
observation, we can expect that the Ath membrane instanton is given by a periodic function. 
Especially, if we assume that it is expressed by using the trigonometric function tan ^ and 
simple rational functions, we can determine the function form of the membrane instantons by 
combining the results of the WKB expansion. In fact, this is how (ll.38j) was first determined 
in [221(35]. 

Note that there are no other exponents e _ ^ +2 ^ p labeled by (m,£) in fll.24p which have 
the same contribution as the first and the second membrane instantons e~ 2/i and e -4p except 
the pure worldsheet instantons. For the higher instanton effects, however, the situation is more 
difficult. The reason is that besides the worldsheet instantons and the membrane instantons 
we need to consider the bound states between them e ~(uf+ 2 ^ ( m ^ 0,£ ^ 0) [23],[35]. We do 
not have clues to them from any systematic expansions because the bound state effects are 
detected neither from the’t Hooft expansions nor from the WKB expansions. The only clue is 
the pole cancellation mechanism as we have explained in section 11.51 and used above to fix the 
coefficients of the membrane instantons. Luckily, we find that if we assume that the bound 
state effects are expressed in terms of the corresponding worldsheet instantons d m (k) 
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(3.16) 


and the corresponding membrane instantons ae(k ) as (A m ^(k) = —2rmi 2 ag{k)) [23] 


(.Pi,- ,Pi) 


<1 

pdp 2 ! • • -pd 


with the sum running over all of allowed partitions (pi, • • • ,Pe) of £ satisfying p 1 + 2 p 2 + • —h 
£pi = £, we can observe a general cancellation among the worldsheet instantons, the membrane 
instantons and their bound states, whose finite results coincide beautifully with the numerical 
studies. This observation implies that, if we redefine the chemical potential p as 

1 OO 

Peff = p +— ^a e (k)e~ 2 ^, (3.17) 

k £=1 

or more explicitly as in (II.23ft for integral fc, 


• the bound state effects are absorbed into the worldsheet instanton effects. 


Besides, we have a bonus from this redefinition: 


the quadratic polynomial coefficients of the membrane instantons are reduced into linear 
polynomials, 


C 


J (peff) = + B ^eS + A+^ d m(k)e + ^(^(/c)/i e ff + c e (k))e 2£/ieff , (3.18) 


m= 1 


£=1 


and the coefficients ce(k) is expressed in terms of be(k) as 

2 d b t (k) 


ce(k) = —k‘ 


dk 2 £k 


(3.19) 


3.5 Topological string theory 

In the above several subsections, we have seen that the WKB expansion, computation of the 
partition function and the pole cancellation mechanism can determine all of the worldsheet 
instantons, the membrane instantons and their bound states order by order. It is not clear, 
however, how the coefficient functions of higher orders generally look like. In this subsection, 
we shall explain a further relation between the ABJ(M) matrix model and topological string 
theory and determine the coefficient functions systematically. 

We shall sketch a relation between the membrane instanton and the refined topological 
string theory in the NS limit after recapitulating the relation between the worldsheet instanton 
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and the (unrefined) topological strings. The large N expansion of the free energy in the 
ABJ(M) theory is formally equivalent to the free energy of pure Chern-Simons theory on KP 3 , 
which in turn is related to the topological string theory on a non-compact Calabi-Yau manifold, 
known as the local P 1 x P 1 [42]. It turns out that the chemical potential in the ABJ(M) Fermi 
gas is a natural variable in the large radius expansion on the topological string side. The 
Kahler parameters of two P^s in local P 1 x P 1 are identified with the chemical potential as in 
(ll.2ip . and they are exchanged under the Seiberg-like duality of the ABJ theory. Surprisingly, 
it was found in [24] that the non-perturbative corrections in the ABJ(M) matrix model are 
completely determined by the refined topological string on the local P 1 x P 1 , and the precise 
relation is given by (1 1 . 16 [) . The first term in (II. 16ft represents the worldsheet instanton, which 
corresponds to the standard topological string. This part of the correspondence (II. 16jl is a 
usual story for the genus expansion of matrix models, and it can be checked, in principle, to 
any order of the string coupling g s = 2 /k. The non-trivial part of the proposal in [24] is the 
second term in (II . 16j) . namely, the correspondence between the membrane instantons in the 
ABJ(M) matrix model and the refined topological string in the NS limit. Note that in (11.16]) 
the pole cancellation mechanism works within every BPS multiplet. Namely, we can check 
that the combination of Atop and Fns in (II . 16j) with common BPS indices A r j] - r does not 
contain any divergences. Also, we stress that not only the singular part but also the finite 
part are correctly reproduced from (11.161) . 

As discussed in |58b60] . the refined topological string theory in the NS limit are closely 
related to the quantization of mirror curves of local Calabi-Yau manifolds. One can easily see 
that the Fermi surface H(x,p ) = 0 of the ABJM Fermi gas corresponds to the mirror curve 
of the diagonal local P 1 x P 1 . More generally, the mirror curve of local P 1 x P 1 is given by 


H(x,p ) = — 1 + e x + e p + zie x + z 2 e p = 0. (3.20) 

We can “quantize” this curve with [x,p] = ih , where the Planck constant h of this quantum 
mechanical system is unrelated to the previous Planck constant in the Fermi gas system and 
is determined later. The corresponding difference equation for the wave function T(a;) is 

(—1 + e x + z\e~ x )^{x) + T(a; + K) + z 2 ^(x — h) — 0, (3.21) 


whose formal solution 

T(a;) = exp 

gives the “quantum periods” 


-S(x,h) 


, S{x 1 h) = Y J Sn{x)h 2 \ 


n =0 


u A = (j) ds, n B = (j) ds. 


(3.22) 


(3.23) 
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Using the quantum periods, we can compute the free energy in the NS limit [60]. 


After identifying the two moduli Zj (/ = 1,2) in (13.2(1 and h as 

Tj 27 TX 

zj = e - ^, h =-= mk , (3.24) 

9s 

we find that the membrane instanton coefficients ae(k) and be(k) appear as the coefficients 
of the quantum A-period II ^ and the quantum B-period II#, respectively, in the small zj 
expansion. Let us see explicitly how it works in the following. 

We shall explain how to compute the quantum periods systematically up to any desired 
order in the small zj expansion. As discussed in [24([60] it is convenient to introduce V(x) by 


V(x) 


ty(x + h) 

'k(x) 


(3.25) 


Then 63.21)1 is rewritten as 


Zl 


V(X) = 1-X-^ 


Z 2 


(3.26) 


X ViXq- 1 )’ 

where we have also introduced X = e x and q = e h . One can easily solve (13.26ft in the small zj 


expansion as 


V(X) = l-X-%+ qZ2 


q 3 z l z 2 


+ 


q A zl 


X X-q X(X- q y (X- q y(X-qi) 


+ 0(z 3 


(3.27) 


Once we know V(X), the quantum periods are obtained by rewriting the contour integral as 
a sum of residues pj 


U A (zt, z 2l q) 

n B (^i,2 2 ;g) 


-Res x=0 — log 


Res 


X=qi ' 


1 - A 
log(X) 


j> o 


X 


log 


v{xy 

l-X 


The quantum A-period is given by 

3 

IIa(^i, z 2 , q) = Z\ + z 2 + (4 + q + q l )z\Z 2 + -(z^ + z%) + 0(z 3 ), 
which defines a quantum version of the mirror map 

llog — = n A (z 1 ,z 2 ,q). 

L Zl 

§Here we will use slightly different normalization of quantum periods as compared to [24] , 


(3.28) 


(3.29) 


(3.30) 
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Here Qj is related to the “effective” Kahler parameter Tf s by 



(3.31) 


which takes care of the effect of the bound states. Inverting the quantum mirror map, we can 
express zi in terms of Qi 



On the other hand, the quantum B-period is related to the derivative of free energy in the NS 
limit as 



where n^ en denotes the even-power part of the quantum B-period in the h expansion. Calcu¬ 


lating the residue in (j3.28j) , the right-hand-side of (13.33P is easily found to be 



(3.34) 


Rewriting the above expression in terms of Qj using the quantum mirror map (I3.32p and 
integrating (13.331) . finally we arrive at 



QiQ 2 + 0(Q 3 ). 

(3.35) 


Fns(Q\, Q 2 , q) 


Plugging this expression into (1 1.16 p . we can show that our formula (ll.lbp reproduces the 
all known results of the membrane instantons and the bound states in the ABJ(M) matrix 
model [1241139] . 

From this non-trivial correspondence in the ABJ(M) theory between the membrane instan¬ 
tons and the refined topological string theory in the NS limit, in [24] it is further conjectured 
that this correspondence holds for general local Calabi-Yau manifolds, and the combination of 
F t op and Fns appearing in (j 1 .16 p gives a non-perturbative completion of the topological string 
theory. However, there is one subtle point: we have to turn on an extra discrete R-field along 
the worldsheet instanton part in order for the pole cancellation to work for the general case. 
See [24] for the detail. 

In |6TM65] . it is argued that this correspondence between the membrane instantons and the 
refined topological string theory in the NS limit in the ABJ(M) theory can be naturally under¬ 
stood from the “dual” perspective, by studying the spectral problem (I2.16P of the Hamiltonian 
in the Fermi gas system. 
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4 Related topics 

In the previous section, we have reviewed the progress in understanding the partition functions 
of the ABJM theory and the ABJ theory. Here we shall explain some progress after finding 
the exact instanton expansion of the ABJ(M) partition function. 

4.1 Exact results with supersymmetry enhancement 

As seen before, the large n expansion of the grand potential is completely determined by the 
refined topological string theory on local P 1 x P 1 [24]. However, it is not technically easy to 
compute the higher order instanton corrections in practice. It is well-known that at k = 1,2, 
the supersymmetry of the ABJM theory is enhanced from J\f = Q to J\f = S. It is natural 
to expect that for these levels some simplifications happen. In fact, as we will see below, the 
large /i expansion of the grand potential drastically simplifies at k — 1,2 [[66j. We can write 
it down in a closed form in terms of the topological string free energy. 

Let us first write down the results in [66]. The complete large /i expansion of Jfc(/i) at 
k = 1,2 is exactly given by 




We need to explain the notations in these equations]^] The functions F 0 (t) and F\ (t) are 


the standard genus zero and genus one free energies on local P 1 x P 1 , respectively. These 
are computed in a standard way of the special geometry. The function F± s (t) is the first 
correction to the refined topological string free energy in the NS limit. The Kahler modulus t 
is related to the complex modulus z by the mirror map 



(4.2) 


The only difference between k = 1 and k = 2 is the identification of z and /r: 


z = e 4/i for k = 1 , z = e 2/i for k = 2. 


(4.3) 


^ We denote the effective Kahler modulus T eS by t in the current case with fixed k. Since there are two 
Kahler moduli t\ and t 2 in local P 1 x f 1 , the free energy is in general a function of these two parameters 
(ti,t 2 ). Here we denote the free energy in the diagonal slice by F g (t,t) = F g (t). 
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At large radius point (t —> oo), the genus zero free energy is expanded as 


f 3 

6 


F 0 (t) = -- 2C(3) + 4e _t - -W 2t + 0(e 


-31\ 


(4.4) 


Eliminating t by (14. 2 j) . one easily finds 

1 


W) - tW) + = -g log 3 o - 2C(3) 

+ 4(log 2 z — logz + 1 )z — ^(52 log 2 z — 2 log z T 9)z 2 + (9(^ 3 ). 


(4.5) 


The free energies E\(f) and F 1 NS (t) are exactly given by 


Fiit) = log[64^(1 + 16z)] - ^ log 


K{—lQz) 


71 


log Z 2 10 n o. 

-k + r~ T 3 +0(z) - 


U NS (i) = ^ >“go - 5 J log ! 1 + 16o) = ^ - fz + fo 2 + 0(z 3 ). 


(4.6) 


where K(z) is the complete elliptic integral of the first kind. Plugging these results into (14.ip . 
one can check that the large /i expansion in table [4] is correctly reproduced. 


In the following, we briefly explain how to obtain these results. Here we focus on the 
case of k — 2. What we should do is to take the limit k —> 2 in the general formula (ll.lOp . 
As mentioned before, both the worldsheet instanton correction and the membrane instanton 
correction are divergent in this limit, but the sum of them is finite. Since we already know the 
divergences are canceled, we can focus on the finite parts. Let us see the worldsheet instanton 
part in fll.lfij) . It is easy to see that the finite part comes only from g — 0,1. The g — 1 part 
is trivial. In the limit k —> 2, each term of the g = 0 part is 


f 2llW \ Ad'Wjj. 

lim 2 sm —-— e k ~ 

k—>2 \ k J 


^—2dwfi 

12n 2 w 2 


(3 + n 2 w 2 + 6 dwg + 6 d 2 w 2 ^ 2 ). 


(4.7) 


where ~ means that we extract only the finite part. Therefore the finite part of the worldsheet 
instanton correction is finally given by 


1 

47r 2 


pr\t) 




inst 


(t) + 


^dtFrv+Frit), 


(4.8) 


where F™ st (t) is the instanton part in F g (t). We have used the relation t = 2 /i e g at k — 2 . 
Similarly, as shown in [ 66 ], the hnite part of the membrane instanton correction is finally given 
by 


1 

8 


(<\ 


dt 2 )' 2 Fo(ti, t 2 )\t 1= t 2 =t + F 1 NS ' lllht (f) 


t 

8 



T-iNS,inst 

w 


( 0 - 


From these results, one obtains the exact expressions (14.11) . 


(4.9) 
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Once the modified grand potential is known, one can compute the grand partition function. 
As seen before, the grand potential is constructed from Jk{p) by (11.81) . with the subscript M 
dropped here. Let us plug the result (14.ip into (11. 8 p . For k — 2, using 


e J 2 (/i+27rm) _ e J 2 (k) eX p 


7T in 2 T + 27 rin 



(4.10) 


with 

2 ? 1 

r = e = ^2 W (t) ~ Km (4-11) 

we find that the grand partition function is expressed in terms of the elliptic theta function 
by 

S2M = e*°W«-l,T). (4.12) 

Similarly, the grand partition function at k — 1 is given by 

S M = e ^ 3 - l, T y (4.13) 


We stress that these expressions are exact. All the instanton corrections are encoded through 
the topological string free energies F 0 (t), F\ (t) and F, NS (f). 

There are two important consequences of the exact expressions (14.12)) and (14.131) . One 
is that one can know the eigenvalues of the density matrix p as the zeros of Eh,(/i). In fact, 
it is well-known that the theta function has an infinite number of zeros. It was confirmed 
in [ 66 ] that the zeros in (14.121) and (14.131) precisely reproduce the eigenvalues of p at k = 1, 2, 
numerically computed by the exact quantization condition in [61] . The other important point 
is that it is possible to analytically continue S^(/i) from p—>-ootop—)■—oo. The continuation 
is essentially the modular S-transformation r — > t — — 1/t. In [ 66 ], this analytic continuation 
was indeed performed, and then the small e M expansion was computed. The results perfectly 
agree with the exact values of the partition function in table [3] (see [2T1I22] for more values)! 

Furthermore, in [66.;,[6Z] the exact results for ABJM at k — 4 and ABJ at (k, M ) = (2,1) 
were presented. These also reproduce the known results perfectly. Note that at k = 4 the 
supersymmetry is no longer enhanced. Nevertheless, we can write down the grand potential 
in a closed form. Moreover, in m, exact functional equations among the grand partition 
functions in the ABJ matrix model were found. These relations are quite similar to so-called 
quantum Wronskia.n relations [ 68 ]. It seems to imply an interesting connection to integrable 
models. 
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4.2 ’t Hooft expansion and Borel resummation 


Throughout this paper, we focus on the large N expansion in the M-theory limit (II.32ft . As 
already seen in subsection 11.41 it is also interesting to consider the ’t Hooft limit (jl.26j) . In 
this limit, the ABJM theory is dual to type HA string theory on AdS^ x CP 3 , and the large 
N expansion is related to the genus expansion in this string theory. In the ’t Hooft limit, 
the free energy admits the large N expansion (ll.27p . The genus g free energy F g ( A) can be 
exactly computed by using the holomorphic anomaly equation order by order G5|. Then, the 
genus expansion (11.271) turned out to be an asymptotic series for given finite A. The standard 
way to resum such formal asymptotic series is known as Borel resummations. For the genus 
expansion (11.271) . let us define its Borel transform by 

mo = £ (2yvyi c2 ^ 2 ' (414) 

Remarkably, as observed in [291169], this Borel transform has no singularities on the positive 
real axis. Therefore the following Borel summation is well-defined: 

r°° A( 

SF(N, A) = N~ 2 F 0 (X) + Fi(X) + / e-tJBF(NC). (4.15) 

Jo C 

The crucial conclusion in [69] is that this Borel summation does not reproduce the exact values 
of F(N, A). There is an exponentially suppressed correction between them: 

F exact (AA A) - SF(N, A) = 0(e~ nN ^). (4.16) 

This correction is nothing but the membrane instanton correction! ft is not captured by the 
perturbative 1/N expansion. Thus in the’t Hooft limit, the perturbative genus expansion and 
its Borel resummation are insufficient to reconstruct the exact free energy. We need to take 
into account the membrane instanton corrections in addition to the perturbative expansion. 
It is an open problem how to incorporate such membrane instantons in the Borel analysis. As 
pointed out in [29] . such corrections should be related to complex singularities of the Borel 
transform (14.141) . 

As reviewed in this paper, we already know the complete M-theoretic expansion of the 
free energy with the help of the topological string. It should be in principle possible to 
use this result to understand the membrane instanton corrections to the genus expansion. 
Interestingly, it was observed in m that the Borel resummation (14.15ft already contains a 
membrane instanton-like correction in the M-theoretic expansion from the Fermi gas result. 
It would be interesting to reveal a quantitative relation between the ’t Hooft expansion and 
the M-theoretic expansion. 
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4.3 Other theories 


In section 14.11 we have seen special cases with supersymmetry enhancements. Here, let us 
see how the analysis can be generalized to other theories of the M2-branes. Especially, when 
we are interested in gaining the universal property of the M2-branes, it is very important to 
generalize to as many theories as possible. We shall list some results obtained recently. 

First, let us give a viewpoint from the quiver diagrams, where the vertices represent the 
gauge group factor U(N) and the edges represent the bi-fundamental matters between the 
two factors connected. Using the quiver diagram, the ABJM matrix model is nothing but the 
Dynkin diagram of affine A\ or A\. Then, it is natural to ask how about a general matrix 
model associated with A r . 

It was already found in [32] that if the density matrix is hermitian, the perturbative part of 
any A r matrix model is summed up to the Airy function with the coefficient C given explicitly 
in [321 ITT] . The other coefficients are not easy to obtain in general. In [00] it was found that 
as long as the quiver diagram associated with the levels is a repetition of an original quiver 
diagram, we can express the grand potential of the new quiver diagram completely in terms 
of that of the original diagram. In [72] it was restricted to the Af = 4 case [73] where the 
Chern-Simons levels are given by 

k 

ka = ~{s a - Sa-i), s a = ±l, (4.17) 

and the explicit form of the coefficient B was found to be 

B (°) 

B=— + hBW, (4.18) 

k 

with 


B {0) 

H (2) 


£(p) | S(g) 
6L £(<?) £(p) 


1 

24 


E(g)E(p) - 12 


4 


£(g)£(p)J’ 

z(q,p,q) 

E(g) 


Z(p,q,p) 


S(g,p)S(p,g) V 

£(g)£(p) ): 


(4.19) 


Here we adopt the notation of E(L), with L denoting an alternating sequence of q and p, 
whose definition is given by 


s (<?) = XA’ Z(p) = J2'Pa, 

a— 1 a =1 

£( q,p)= qaPb ’ s (^’^) 

l<a<b<m 


l<a<b<m 
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PaqbPc- 


(4.20) 


£(g,p,g)= QMc Y>(p,q,p)= 

l<a<b<c<m 1 <a<b<c<m 

Note that the condition in each sum can be stated as the requirement that we choose q a and 
p a out of a sequence qi,pi, qiiPii ■ • • , qm,Pm by respecting its ordering. 


Furthermore, if we restrict ourselves to the special AT = 4 (q,p) model where the q edges 
with s a = +1 and the p edges with s a = — 1 are separated, we find that the coefficient A is 
given in terms of that of the ABJM matrix model by [72] 

A = ~(p 2 A\Bm(qk) + q 2 Aabjm(p^)) • (4-21) 

A similar proposal was made in [4T] for the so-called M-theoretic matrix model [551174] , whose 
Fermi gas Hamiltonian is related to that of the (q, 1) model by a canonical transformation [75]. 
In [76] . the membrane instanton of the (q,p) model was studied and a systematic expansion 
in terms of the derivative operators acting on the classical expression was found. Using the 
results we were able to study the ( 2 , 2 ) model in full details [[77] and find that the non- 
perturbative effects are given by the refined topological string theory as in the ABJM matrix 
model (11.121) . It still awaits to identify the Calabi-Yau manifold whose BPS index describes 
the ( 2 , 2 ) model, but we observed that the diagonal Gopakumar-Vafa invariant is nothing but 
that of the ZJ 5 del Pezzo surface. It is interesting to find that the correspondence between the 
topological string theory and the matrix model coming from the superconformal theory, is not 
restricted to the case of the ABJM theory. In [78], the non-perturbative effects were further 
explored. Based on the systematic WKB analysis and exact computations in various special 
cases, analytic expressions of the first few membrane instanton corrections and worldsheet 
instanton corrections were conjectured for general (q,p) and k. 

The generalization, however, is not restricted to the A quiver. In [79] an ADE sector 
was specified by requiring the long range force among the eigenvalues vanishes. It is then 
interesting to ask whether there is a Fermi gas formalism for these series and whether the 
perturbative corrections are still summed up to the Airy function. In [SDUSI] we continue to 
study the D quiver. We find that after using the Cauchy determinant (12. 2\i in a different 
way we can still give a Fermi gas formalism to this type of theories and prove the previous 
proposal [82] 


C 


7 r 2 k \cr 0 CTi 


r—1 

£ 

m—1 


Sm ^m+1 
&m&m +1 



(T r (T r _l_l J 


(4.22) 


with the variables a m given as 


r r 

^ |Sm &n\ \^m ^n|) 4|s m |, (Jq 2 (t 2), ^Y+l 2 ^ ^ l^nl? (4.23) 

71=1 77=1 
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by the reordered s, 


0 < |s r | < |s r -i| < ••• < |si|. (4.24) 

As in the Af = 4 case of the A quiver, here we can also restrict ourselves to the case of uniform 
s a , s a = 1, and find that the coefficient A is consistent with 

A = -(AABjM(2r/c) + r 2 AABjM(2(r — 2 )fc)), (4-25) 

up to 0(k 5 ), which is very similar to (14.21)1 for the (q,p) model. We hope that all these 
findings will finally be completed to a more systematic understanding. 


5 Summary and further directions 

In this review article, we have explained carefully the quantum effects of the ABJM theory 
and the ABJ theory, where the partition function and the vacuum expectation values of the 
half-BPS Wilson loop can be studied in full details including both the perturbative and non- 
perturbative corrections. 

We have seen many interesting properties. For example, the leading large N behavior of 

1 3 

the free energy is given by ~ k^Nz. The perturbative 1/N corrections are summed up to an 
Airy function. The non-perturbative corrections include the effects interpreted as worldsheet 
instantons and those interpreted as membrane instantons. In general, the coefficients of both 
instantons and their bound states can be divergent at certain values of k, though as a whole 
the divergences in these instanton effects cancel among themselves and the final result is 
completely finite. 

Also, in section 14.31 we have seen that many generalizations enjoy these interesting proper¬ 
ties. Hence, one further direction is to generalize our analysis to as many theories as possible 
and identify this class clearly. It is surprising to find that the superconformal Chern-Simons 
theory of the D quiver can also be formulated with the Fermi gas system. It is then natural 
to ask whether we can apply the same method even for the E quiver or their cousins with the 
orthogonal groups or the symplectic groups [HlSIj. We hope that through these studies we 
will end up with a universal understanding of the M2-branes. 

Next, let us also stress that, after a large number of numerical analysis, our large /i expan¬ 
sion of the grand potential (or large N expansion of the partition function) remains to be a 
conjecture. It is of course interesting to see how this conjecture is proved from a mathematical 
method. The study explained in section |4~T1 would be a first attempt in this direction. 
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Finally, let us come back to a more physical question. Our original motivation in studying 
the AB JM partition function is to understand the M2-branes. After seen all of the perturbative 
and non-perturbative corrections, it is reasonable to come back to ask what we have learned on 
the M2-brane interactions. Unfortunately, we do not have much to say about it at this stage, 
except that some of the results can be rederived from the supergravity [MHSB] . However, we 
believe that the beautiful structure we find would be a guide towards a clearer understanding 

3 

of the Nz law or even a better description of the mysterious M-theory. 
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